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Abstract 

We discuss a novel kind of nonlinear coupler with one channel filled with a negative index 
i material (NIM). The opposite directionality of the phase velocity and the energy flow in the NIM 

channel facilitates an effective feedback mechanism that leads to optical bistability and gap soliton 
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Nonlinear optical couplers have attracted significant attention owing to their strong po- 
tential for all-optical processing applications, including switching and power-limiting devices. 
Transmission properties of a nonlinear coherent directional coupler were originally studied 
by Jensen [1] , who concluded that a coupler consisting of two channels made of conventional 
homogeneous nonlinear materials is not bistable. 

Bi- (or multi-) stability is a phenomenon in which the system exhibits two (or more) steady 
transmission states for the same input intensity [21 E]- Optical bistability has been predicted 
and experimentally realized in various settings including a Fabry-Perot resonator filled with 
a nonlinear material [3] and layered periodic structures [I] . In this Letter we describe a novel 
nonlinear optical coupler structure that utilizes a negative index metamaterial (NIM) p2, E] 
in one of the channels and a conventional positive index material (PIM) in another channel 
as shown in Fig. [TJ The linear transmission properties of a similar optical structure were 
previously studied by Alu and Engheta [7]. We show that such nonlinear coupler (NLC) 
can be bistable. Bistability occurs owing to the effective feedback mechanism enabled by a 
fundamental property of NIMs - opposite directionality of the wave vector and the Poynting 
vector. Moreover, our results suggest that the entirely uniform PIM-NIM coupler structure 
supports gap solitons - a feature commonly associated with periodic structures [I], [8]-|14j. 
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FIG. 1: A schematic of a nonlinear PIM-NIM coupler. Light is initially launched into channel 1 
(PIM). A wave vector and a Poynting vector are parallel in the PIM channel and anti-parallel in 
the NIM channel enabling a new backward-coupling mechanism. 



Continuous wave propagation in a nonlinear coupler can be described by the following 
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system of equations: 
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where a\ and a 2 are the complex amplitudes of the modes in the PIM and NIM channels 
respectively, K\ 2 and n 2 \ are the coupling coefficients defined as in Ref. [T|, S = Pi — f3 2 
is the mismatch between the propagation constants in the individual channels, /3| = 
(lu / c) 2 fij(uJo)e j(^ ) , (j = 1, 2), Oj is the sign of the refractive index rij = a/ ^^ujqis^ujq), Ej 
and /ij are linear frequency-dependent dielectric permittivity and magnetic permeability, u 
is carrier frequency, x) is the nonlinear (electric) susceptibility, and c is the speed of light 
in a vacuum. In the case of PIM-NIM coupler o\ is positive, while o 2 is negative. 

Assuming the form ai j2 = U\ j2 exp(iqz) exp(^fiSz/2) for the solutions of Eqs. ([TJ [2| in the 
linear regime, we find the following relation between q and 5 for the PIM-NIM coupler 

q 2 = (S/2) 2 - k 12 k 21 , (3) 

which indicates the presence of a bandgap for \5\ < 2 v /k 12 /«2i- The photonic bandgap 
is usually a feature typical for periodic or distributed feedback (DFB) structures such as 
fiber Bragg gratings (FBGs) or thin film stacks [7j. Formation of the bandgap in a uniform 
structure considered here is one of the unique properties of the PIM-NIM coupler arising from 
introduction of the NIM into the nonlinear coupler. Introducing two parameters a 2 = uf + u 2 . 
and / = u 2 /ui, the nonlinear counterpart of the relation [3] can be written in the form 

n 21 + f 2 K 12 a 2 (7i +72/ 2 ) 
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where 
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is the nonlinearity coefficient. Figure [2] shows linear (a) and nonlinear 5 — q curves (b)-(d) 
for the case of self-focusing Kerr nonlinearity. In Fig. [2^b) both PIM and NIM channels are 
nonlinear, with the same nonlinearity coefficients. Beyond a critical power level, the lower 
branch of the 5 — q curve forms a loop. Both linear and nonlinear 5 — q curves in Fig. |2^a) 



and[2](b) resemble dispersion relations found in the case of linear and nonlinear FBG [8| [15] . 
with the difference that in the nonlinear case the effect of cross-phase modulation has been 
neglected and forward and backward waves are spatially separated. However, importantly, 
in the case of PIM-NIM NLC both channels are made of homogeneous material with no 
periodicity or an external feedback mechanism (such as a cavity). However, the effective 
feedback mechanism is provided by the inherent property of the NIMs, that is opposite 
directionality of the phase and energy velocities. As shown in Fig. [TJ while the propagation 
vectors of the waves propagating in both NLC channels point in the same direction assuring 
a necessary phase-matching condition, the Poynting vectors corresponding to the energy 
flow direction point in opposite directions. As light propagates in the PIM channel in the 
forward direction it continuously couples to the NIM channel where it flows in the backward 
direction. Therefore, the PIM-NIM coupler acts as an effective DFB structure. Figure |2](c) 



corresponds to the case of different nonlinear coefficients in the two channels, while Fig. \2, 



corresponds to the case of the nonlinear PIM and linear NIM channel. Figures p[c) andg 



also illustrate that there are more degrees of design freedom in the NLC case in comparison 
with the DFB structures since in general case the nonlinear coefficients 71 and 72 as well as 
the coupling coefficients K12 and K21 may be not identical and can be properly designed. 

Although in general case K12 7^ ^21 and 71 7^ 72, in order to illustrate basic new physical 
effects associated with the PIM-NIM NLC, in the following discussion we assume identical 
linear coupling coefficients K12 = K21 = K and nonlinear coefficients 71 = 72 = 7. Then, 
making the substitution a± = A 1 exp(z0 1 ) and a 2 = A 2 exp(i(j) 2 ), where A X) A 2 , 4>i and </> 2 
are real functions of z, the equations for A\, A 2 and ip = (pi — (f) 2 + 5z can be written in the 
form 
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dz 




From the system of equations (pj) the constants of the motion are given by 
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FIG. 2: PIM-NIM coupler 5 — q relations: (a) linear (k\i = Ki\ = 5), (b) two channels with the 
same nonlinear susceptibility 71 = 72 = 7 (K12 = K 2i = 5, a 2 7 / 'k = 6), (c) two channels with 
different nonlinear susceptibilities (kyi = k>%\ = 5, « 2 7i/k = 6, oP'^/k = %)> (d) the NIM channel 
is linear, while the PIM channel is nonlinear (k% 2 = K21 = 5, a 2 7i/ 'k = 6, a 2 ^/^ = 0). 

where Pi = A\ and P2 = A\. C is defined by the boundary conditions at z = and z = L, 
where L is the length of the coupler. The expressions ^ should be compared to that of 
conventional PIM-PIM NLC, in which case C = P x + Pa- 
Then, power evolution in channel 1 is described by the equation 



{dPt/dzy = k(4k + r 7 )P!P 2 + Tk5P 1 - k 2 T 2 /A -(5 + 7P2) ^1 • 



(9) 



If light initially is launched into channel 1, i.e. Ai(0) = A , A 2 (L) = 0, then C = A 2 (L). 
In the case of 5 = and T = Eq. ^ reduces to 

(dPjdz) 2 = Ak 2 P 1 (P 1 - C) - 7 2 P! 2 (Pi - C) 2 . (10) 



Introducing new variable Z = Cuaz and y 2 = a 2 P\(Pi — C), where a = j/2k, Eq. (10) can 
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be written in the form 

(dy/dZ) 2 = (l-y 2 )(l + k 2 y 2 ), (11) 



where k = (aC/2) 1 . The solution of Eq. (11) is given by 



y(Z)=m{Z-Z 1 ;ik), (12) 
where Z\ = CuaL as follows from the boundary condition, sn(z',k') is the Jacobi elliptic 



function [16] . Using the properties of Jacobi elliptic functions and Eq. ( 12 ) the solutions for 
Pi and P2 can be written in the form 

_ ^ l + dn[2/<z-L)/m;m] 
Fl{Z) ~ ° 2dn[2/<* -L)/m;m] ' Udj 
PM = r l-dn[2K(z-L)/m;m] 
21 J 2dn[2K(z-L)/m;m] ' 1 J 

where 

fc 1 

m 



v / TTP ^/1 + ( 7 C/4k) 2 ' 
The parameter C can be found using the transcendental equation 



, 1 + dn[2/tL/m; m] 



2dn[2 K L/m;m] ' (15) 
Finally, one can define transmission and "reflection" coefficients for the nonlinear coupler 



as 



/ML) = C^ 2dn[2/ t L/m;m] 

Pi(0) A 2 1 + dn[2«L/m; m] ' 1 J 

«t 1 - *W± = 1 ~ dn[2/tL/m; m] 

Pi(0) l + dn[2«i/m;m]' 1 J 

Figure |3](a) shows output power Pi(L) as a function of input power -Pi(O) for three values 
of kL = 2 (solid line), kL = 4 (dot-dashed line), and kL = 6 (dashed line), assuming k is 
varying, the coupler length L = 1 (in the units of length) is fixed and 7L = 6. As the coupling 
between the channels increases and the effective feedback mechanism establishes, PIM-NIM 
NLC becomes bistable or more generally multistable such as in the case of kL = 6. Its 
transmission characteristics are very similar to those of the DFB structures [E]-[lI] with an 
important fundamental difference that bistability in the PIM-NIM coupler is facilitated by 
the effective feedback mechanism originating from the NIM's intrinsic property. Figure ^b) 
shows the transmission coefficient 1 as a function of input power -Pi(O) for two values of 
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jL = 1 (solid line) and 7L = 6 (dashed line), assuming 7 is varying, the coupler length L = 1 
is fixed and kL = 6. As the nonlinearity coefficient decreases the threshold of bistability 
shifts to the higher values as expected. 




Input Power PrfO) Input Power P^O) 

FIG. 3: (a) Output power Pi{L) from channel 2 as a function of input power Pi(0) launched into 
channel 1 for three values of kL = 2 (solid line), kL = 4 (dashed line) and kL = 6 (dot-dashed 
line) when 7L = 6. (b) Transmission coefficient defined as P\{L)/P\{Q) as a function of -Pi(O) for 
^L = 1 (solid line) and 7L = 6 (dashed line) when kL = 6. Transmission resonances are indicated 
by the numbers 1,2,3, and 4- 

The phenomenon of bistability in DFBs is closely related with the notion of gap soli- 
tons [7]- [15]. As shown in Fig. |3](b), the transmission coefficient approaches % = 1 at 
the points 1, 2, 3 and 4, suggesting the existence of transmission resonances. At the reso- 
nance, corresponding to a point 1 in Fig. [3|b), spatial power distributions P\{z) (solid line) 
and P2{z) (dashed line) peak in the middle of the structure as shown in Fig. |lj The dot- 
dashed line in Fig. [4] shows the constant of the motion C = P\ — P2. At this transmission 
resonance incident light is coupled to a soliton-like static entity that peaks in the middle of 
the structure and is known as a gap soliton pTJ| [TTj . Nonstationary gap solitons in PIM-NIM 
NLC analogous to those found by Aceves and Wabnitz [T3] in a context of periodic media 
will be discussed elsewhere. It is notable that a gap soliton, usually existing in periodic 
structures, forms in a uniform structure in the case of the PIM-NIM coupler, owing to the 
new backward-coupling mechanism in the PIM-NIM coupler. In summary, we found that 
backward coupling between the modes propagating in the PIM and NIM channels enabled 
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Distance 

FIG. 4: Spatial power distribution in channel 1 (solid line), in channel 2 (dashed line), and the 
constant of motion C = P\ — P2 (dot-dashed line) versus z at transmission resonance indicated by 
number 'V in Fig. \^b). 

by the basic property of NIMs, oppositely directed phase and energy velocities, results in op- 
tical bistability in PIM-NIM NLC and gap soliton formation. These effects have no analogies 
in conventional PIM-PIM couplers composed of uniform (homogeneous) waveguides with no 
feedback mechanism. 

This work was supported in part by the ARO Award # 52730-PH, ARO-MURI Award 
# 50342-PH-MUR, by the NSF Award # DMS-050989, by State of Arizona grant TRIF 
(Proposition 301) and by the RFBR Award # 06-02-16406 



[1] S. M. Jensen, IEEE J. Quantum Electron. QE-18, 1580 (1982). 

[2] C. Bowden, M. M. Ciftan, and H. Robl, Eds. Optical Bistability, (Plenum, New York, 1981). 
[3] H. M. Gibbs, Optical Bistability, (Academic Press, Orlando, FL, 1985). 



8 



[4] H. G. Winful, J. H. Marburger, and E. Garmire, Appl. Phys. Lett. 35, 379 (1979). 
[5] V. G. Veselago, Soviet Physics Uspekhi, 10, 509 (1968). 

[6] V. M. Shalaev, V. M., Optical Negative-Index Metamaterials, Nature Photonics 1, 41 (2007). 

[7] A. Alu and N. Engheta, in Negative-Refraction Metamaterials ed. G. V. Eleftheriades and K. 
G. Balmain (Wiley, New York, 2005). 

[8] C.M. de Sterke and J.E. Sipe, in Progress in Optics, edited by E. Wolf (North-Holland, Ams- 
terdam, 1994), Vol. XXXIII, pp. 203-260 

[9] H. G. Winful, Ph.D. Dissertation, University of Southern California, 1981. 
[10] W. Chen and D. L. Mills, Phys. Rev. Lett. 58, 160 (1987). 
[11] D. L. Mills and S. E. Trullinger, Phys. Rev. B 36, 6269 (1987). 
[12] D. N. Christodoulides and R. I. Joseph, Phys. Rev. Lett. 62, 1746 (1989). 
[13] A. B. Aceves and S. Wabnitz, Phys. Lett. A 141, 37 (1989). 

[14] B. J. Eggleton, R. E. Slusher, C. M. de Sterke, P. A. Krug, and J. E. Sipe, Phys. Rev. Lett. 
76, 1627 (1996). 

[15] G. P. Agrawal, Applications of Nonlinear Fiber Optics (Academic Press, San Diego, 2001). 
[16] M. Abramowitz and I. A. Stegun (Eds.), Handbook of Mathematical Functions with Formulas, 
Graphs, and Mathematical Tables, (9th printing. New York: Dover, 1972). 



9 



